arXiv:l 504.05705V 1 [cs.NA] 22 Apr 2015 


CONVERGENCE OF A FINITE DIFFERENCE SCHEME TO 
WEAK SOLUTIONS OF THE SYSTEM OF PARTIAL 
DIFFERENTIAL EQUATION ARISING IN MEAN FIELD GAMES 

YVES ACHDOU * AND ALESSIO PORRETTA t 

Abstract. Mean field type models describing the limiting behavior of stochastic differential 
games as the number of players tends to +cxd, have been recently introduced by J-M. Lasry and 
P-L. Lions. Under suitable assumptions, they lead to a system of two coupled partial differential 
equations, a forward Bellman equation and a backward Fokker-Planck equations. Finite difference 
schemes for the approximation of such systems have been proposed in previous works. Here, we 
prove the convergence of these schemes towards a weak solution of the system of partial differential 
equations. 
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1. Introduction. Mean field type models describing the asymptotic behavior 
of stochastic differential games (Nash equilibria) as the number of players tends 
to +00 have recently been introduced by J-M. Lasry and P-L. Lions [HISllll], 
and termed mean Geld games by the same authors. Related ideas have been de¬ 
veloped independently in the engineering literature by Huang-Caines-Malhame, see 
for example [20) . For brevity, the acronym MFG will sometimes be used for mean 
Geld games. Examples of MFG models with applications in economics and social 
sciences are proposed in mn]. 

The simplest MFG model lead to systems of evolutive partial differential equa¬ 
tions involving two unknown scalar functions: the density of the agents in a given 
state X G namely m = m{t,x) and the value function u = u{t,x). The present 
work is devoted to finite difference schemes for the systems of partial differential 
equations. Although the methods and the theoretical results obtained below can 
be easily generalized, the present work focuses on the two-dimensional case for the 
following reasons: 1) the one dimensional case is easier and allows too special ar¬ 
guments; 2) in dimension two, the description of the discrete methods discussed 
below remain fairly simple. Besides, several important applications of the mean 
field games theory are two-dimensional, in particular those related to crowd dy¬ 
namics. 

In the state-periodic setting, typical MFG model comprises the following system 
of partial differential equations in (0, T) x 

dll 

— (t, x) — i^Auit, x) + H{x, Vu{t, x)) = F{m{t, x)), (1-1) 

(x,Vn(t,x))) =0, (1.2) 

with the initial and terminal conditions 

u{0,x) = uo{x), m(T,x) = mT{x), in T^, (1-3) 

given a cost function uq and a probability density my- 

Here, we denote by = [0,1]^ the 2—dimensional unit torus, and by A, V and 
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dm 


{t,x) + vAmit, x) + div m{t,x) 


dt 


dH 


dp 
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div, respectively, the Laplace, the gradient and the divergence operator acting on 
the state variable x. The parameter v is the diffusion coefficient. Hereafter, we will 
always assume that > 0. The system also involves the scalar Hamiltonian H{x,p), 
which is assumed to be continuous, convex and regular with respect to p. The 
notation ^{x,q) is used for the gradient of p i—?> H{x^p) at p = q. Finally, in the 
term F{m{t,x)), F is a continuous real valued function defined on K+. Hereafter 
the notation Q will be used for the space-time cylinder (0,T) x T^. 

We have chosen to focus on the case when the cost uq depends directly on x. 
In some realistic situations, the final cost may depend on the density of the players, 
i.e. u\t=o = d>o[w|t=o](3;), where $o is an operator acting on probability densities, 
which may be local or not. We will not tackle this aspect, in order to keep the 
discussion as simple as possible. Similarly, by working on the torus T^, we avoid 
the discussion of the boundary conditions, but other boundary value problems with 
for example Dirichlet or Neumann conditions could be considered. It is also pos¬ 
sible to consider different initial conditions than in (IE3I: if there is a condition of 
the type m{t = 0, •) = mg instead of u{t = 0,-) = uq, then the system models a 
planning problem, see |25) for a description of the model and mathematical results, 
and nanz] for new existence and uniqueness results. 

System (ICT-lfT^ consists then of a forward Bellman equation coupled with a 
backward Fokker-Planck equation. The forward-backward structure is an impor¬ 
tant feature of this system, which makes it necessary to design new strategies for 
its mathematical analysis (see |231124) 1 and for numerical approximation. The main 
results on the mathematical analysis of (in])-(lL2l) are contained in the pioneering 
articles [23l[24], but many important aspects of the theory developed by J-M. Lasry 
and P-L. Lions on MFG are not published in journals or books. They can never¬ 
theless be found in the videos of the lectures of P-L. Lions (in French) at College 
de France: see [25]. A very good introduction is also given in the notes by P. 
Cardaliaguet, |5], with a special emphasis on the deterministic case, i.e. = 0 in 
(fTTD-ira. The survey of Gomes et al (T^j also addresses interesting extensions of 
the model, and the so-called master equation first introduced in [25) . 

Depending on the data and on F and H, different notions of solutions can 
be relevant for (fTTIl-lfOl): indeed, if the right hand side of dm is replaced by 
$[m(t, •)](x) where $ is a nonlocal smoothing operator, mapping probability mea¬ 
sures on to functions, if FI depends smoothly on x and if the data uq and 
ttzt are smooth, then classical solutions can be found, see |23l[24]. The same is true 
if e.g. F[ is Lipschitz continuous w.r.t. its second argument p and F in dm is a 
continuous function. The situation is different in the case when FI has a strictly su- 
perlinear growth with respect to p and F is a continuous function: in this case, one 
has to look for weak solutions, see [23] and the recent article |5S] which is devoted 
to weak solutions to Fokker-Planck equations and to the system (fLTt-lfOl). 

Since the (semi-)analytic solutions of the MFG system do not exist in general, any 
attempt to apply MFG models and to get qualitative/quantitative information from 
them must rely on numerical simulations and scientific computing. Therefore, the 
research has also been active on numerical methods for approximating (inD-(IT3]): 
a numerical method based on the reformulation of the model as an optimal control 
problem for the Fokker-Planck equation with an application in economics was pro¬ 
posed in m- Discrete time, finite state space mean field games were discussed in 
m- We also refer to [nnn] for a specific constructive approach when the Hamil¬ 
tonian is quadratic. Finally, semi-Lagrangian approximations have been studied in 

[ndQ]. 

The finite difference method described and studied below has first been proposed 
and discussed in mm- It will be reviewed in § [2] The numerical scheme basically 
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relies on monotone approximations of the Hamiltonian and on a suitable weak for¬ 
mulation of the Fokker-Planck equation. It has several important features: 

• existence and uniqueness for the discretized problems can be obtained by sim¬ 
ilar arguments as those used in the continuous case, 

• they are robust when z/ —>■ 0 (the deterministic limit of the models), 

• bounds on the solutions, which are uniform in the grid step, can be proved 
under reasonable assumptions on the data. 

A first result on the convergence to classical solutions was given in [^. The issue 
of convergence was studied with more details in HDP: in these works, the starting 
point/assumption was the existence of a classical solution of (If.Ill - (If.31) . The proof 
of convergence mainly consisted in plugging the classical solution into the system 
of equations arising from the finite difference method, and use the consistency and 
stability properties of the scheme in order to get estimates and pass to the limit. 

In the present work, the goal is different: we wish to prove that as the grid steps tend 
to zero, the solution of the discretized MFG system converges to a weak solution of 
(mil-do]), without assuming the existence of the latter; so this work will supply as 
a by-product a new strategy for proving the existence of weak solutions. One key 
step will be to obtain a priori estimates on the solutions to the discrete systems, 
and these will mainly come from the fact that the structure of the MFG system 
is preserved by the chosen finite difference method. This step will be achieved in 
§ m Note that § 14.31 is concerned with a priori estimates for the discrete version 
of the Fokker-Planck equation dm. These estimates may have their own interest, 
independently from MFG models, and may be put in relation with recent works 
of Gallouet et al m in the context of finite volume methods. Once these esti¬ 
mates are obtained, the most important difficulty will be to pass to the limit in 
the discrete Bellman equation. The strategy to that purpose is to first prove some 

compactness of the sequence of state-gradients of the discrete solutions, then to 
adapt some techniques proposed by Boccardo, Murat and Fuel, see [8], for studying 
weak solutions of dni)- This will done in S 15.11 and l6l Here also, we think that the 
passage to the limit in the discrete Bellman equation may have an interest for itself. 

2. Finite difference schemes. In the present paragraph, we discuss the finite 
difference method originally proposed in [S]. 

Let Nt be a positive integer and At = T/Nt, tn = nAt, n = Let 

Tl be a uniform grid on the torus with mesh step h, (assuming that \/h is an 
integer Nh), and Xij denote a generic point in T^. The values of u and m at 
are respectively approximated by and mfj. Let m" (resp. m") be 
the vector containing the values u'^j (resp. mfj), for 0 < i,j < Nh indexed in 
the lexicographic order. Hereafter, such vectors will be termed grid functions on 

or simply grid functions. For all grid functions z, all i and j, we agree that 

= ^(i mod Nh),(j mod n^)- 

Elementary finite difference operators. Let us introduce the elementary finite 
difference operators 




and {D 2 u)ij = 


( 2 . 1 ) 


and define DhU as the grid function with values in 1 

{Dhu)i,j = (LI^w)i,j 


( 2 . 2 ) 


Let be the collection of the four possible one sided finite differences at Xij: 

[VftMjiy = u)ij, {D^u)i-ij, {D^u)ij, {D^u)ij-i^ e R^. (2.3) 




We will also need the standard five point discrete Laplace operator 




^2 ij i)' 


For a set ?; = {v'^)n=o,...,NTJ where is grid functions on T^, it will be convenient 
to define the family of grid functions: 


dt.AtV 


yn+l_yT 

At 


n=0,..., Nx — 1 


(2.4) 


Numerical Hamiltonian. In order to approximate the term H{x,'\/u) in (ll.l[) . 
we consider a numerical Hamiltonian g : —>■]&, (x, gi, <72,93,94) H> 9 (x, 91 ,92, 93,94)- 

Hereafter we will often assume that the following conditions hold: 

(gi) monotonicity: g is nonincreasing with respect to 91 and 93 and nondecreasing 
with respect to 92 and 94. 

(g2) consistency: 9 (x, 91,91,92,92) = Lf(a;, 9), Vx G T^, V9 = (91,92) G 
(gs) regularity: g is continuous and of class w.r.t. (91,92,93,94). 

(g 4 ) convexity : (91,92,93,94) g (x, 91 , 92,93,94) is convex. 

We will approximate H{-,'Vu){xij) by g{xij, [^hu]ij). 

Standard examples of numerical Hamiltonians fulfilling these requirements are pro¬ 
vided by Lax-Friedrichs or upwind schemes, see [S]. For Hamiltonians of the form 
H{x,p) ='H{x) + \p\^, j 3 G ( 1 , 00 ), we may choose 

g{x,q) = n{x) + G(9f,9^,93',9|), (2.5) 


where, for a real number r, = max(r, 0) and r = max(—r, 0) and where G : 
(]R+)^ —M+ is given by 


G{p) = \pf = {pI+pI+pI+pI)''^. (2.6) 


Discrete Bellman equation. The discrete version of the Bellman equation is 
obtained by applying a semi-implicit Euler scheme to (ED, 

_ yn 

+ g{xij, (2.7) 

for all points in and all n, 0 < n < TVy, where all the discrete operators have been 
introduced above. Given (m")„=o,...,AfT-i, (ED and the initial condition uF = 
uoixij) for all {i,j) completely characterizes (u")o<n<AfT- 

Discrete Fokker-Planck equation. In order to approximate equation ED , it is 
convenient to consider its weak formulation which involves in particular the term 



m-^[x. vu) 
op 


w{x) dx. 


By periodicity. 



dH, ^ , 

m-^(x. VU) 

op 


w{x) dx 


f dH 

/ m(x)-^(x, Vu(x)) • Vw(x) dx 
Jt2 Op 


holds for any test function w. The right hand side in the identity above will be 
approximated by 


-h- 




9(xjj, [V/iu]ij) • [Vhw]ij = h'^y^Ti,j{u,m)wtj 


^^3 


^^3 
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where the transport operator T is defined as follows: 


^ 


V 


dg dg 


The discrete version of equation is chosen as follows: 


771^"^^ — TTi^ ■ 

^ + iy{Ahm'^)ij +7Ij(M"+\m") = 0, 


\ 

( 2 . 8 ) 

/ 


(2.9) 


for all n = 0,... ,Nt — 1- This scheme is implicit w.r.t. to m and explicit w.r.t. 
u because the considered Fokker-Planck equation is backward. Given u this is a 
system of linear equations for m. We introduce the compact and convex set 


K-h — {(^2 J )o<i,J<A:';i • ^ ^ ^ ^i,j — (^*^^) 


which can be viewed as the set of the discrete probability measures. It is easy to 
see that if m" satisfies (12.91) for 0 < n < Nt and if S ICh, then m” G ICh for 
all n, 0 < n < Nt- 

Remark 1. An important property ofT is that the operator m i—> [—v{Ahm)i^j — 
Ti,j{u,m)) ^ ^ is the adjoint of the linearized version of the operator u i—>■ [—i'{Ahu)ij + 

g{xi,j, [Nhu\ij)'j ^ y 

This property implies that the strueture of ir7)) -fOI) is preserved in the discrete 
version CT-([ID- In partieular, it implies the uniqueness result stated in Theo¬ 
rem hi. ill below. 

Summary. The fully discrete scheme for system (ini),(IIlD,(lI31) is therefore the 
following: for all 0 < i, j < Nh and 0 < A: < Nt 


- v(Af^u^+'^),^j +g{x^,J, 

^ + v{Ahm^)ij +7)j(M''+7m'=) 


At 

with the initial and terminal conditions 

1 r 


= uo{x^,J), 


Nx 

- fj2 


\x—Xij I oo'^h/2 


mT{x)dx, 


= Pimlj), 

= 0 , 


( 2 . 11 ) 


0<t,j<A^;,. (2.12) 


The following theorem was proved in (using essentially Brouwer’s fixed point 
theorem and estimates on the solutions of the discrete Bellman equation): 

Theorem 2.1. Assume that hold, that uo is continuous on and 

that mT G L^(T^) is a probability density, i.e. mr > 0 and fj .2 mT(x)dx = 1; then 
i2.11\} - l2.1A) has a solution such that m" G K-h, Vn. 

Since (I2.11I) - (I2.12I) has exactly the same structure as the continuous problem 
(frT])-(fT3]). uniqueness has been obtained in with the same arguments as in [23] : 

Theorem 2.2. Assume that (g,i)-(s 4 ,) hold and that F is nondecreasing then 
\2.11\) - i2.12\} has a unique solution. 

Remark 2. Efficient algorithms for solving system i2.11\) - i2.12\) require spe¬ 
cial efforts, essentially because of the forward-backward structure already discussed 
above. We refer to m for the description of possible algorithms and numerical 
results. 
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3. Running assumptions and statement of the main resnlt. We now 

summarize the assumptions that will be made in the whole work. 

• uq is a continuous function on 

• rriT is a nonnegative function in such that fj 2 mT(x)dx = 1 

• F is a continuous function on which is bounded from below. 

• The Hamiltoniaifl (x,p) i—!► H{x,p) is assumed to be convex with respect 
to p and regular w.r.t. x and p. 

• The discrete Hamiltonian g satisfies (gi)-(g 4 ) and the further assumption 

(gs) There exist positive constants ci, 02 , 03,04 such that 

9 q{x,q) -q- g{x,q) > Ci\gg{x,q)\'^ - C 2 , (3.1) 

\gq{x,q)\<C3\q\+Ci. (3.2) 


Take for example g as in (12.51) (12.61) . It is clear that gq{x, q)-q = l3G{q-^ > 9 ^ > 93 > 94 ")) 
h.eace gq{x,q) ■ q - g{x,q) = {P - l)G{q^,q^,q^,qt) - Hix). Since \gq{x,q)\'^ = 

/3^ {G{qi ,q 2 ,93 ,qt)Y^^ ^ we see that (gs) holds if 1 < /3 < 2 . 

We can now state the main result of this article, which establishes the conver¬ 
gence of the solutions of the finite difference scheme towards a weak solution of the 
continuous mean field games system. 

Theorem 3.1. Let ivL), (m") he a solution of the discrete system i2.11\) - i2.1S\) 
and Uh,At, ixih^At be the piecewise constant functions which take the values and 
mfj, respectively, in {tn,tn+i) x {ih — h/2,ih + h/2) x {jh — h/2,jh + h/2). There 
exists a subsequence of h and At (not relabeled) and functions u, rfi, which belong 
to L“(0, T; W^’“(T^)) for any a £ [!> |)j such that Uh^At u and uih^At rh in 
L^{Q) for all [3 € [1,2), and {u, fh) is a weak solution to the system in 

the following sense: 

(i) H{-,Du) £ L^{Q), rhF{rh) £ L^{Q), fh[F[p{-,Du)-Du — F[{-,Du)] £ L^(Q) 

(a) {u,fn) satisfies in the sense of distributions 

(Hi) u,fh£ (^^([O, T]; L^(T^)) and u\t=o = uq, rh\t=T = wt • 

Remark 3. We recall that, if F is nondecreasing and p 1 — F[(x,p) is strictly 
convex at infinity, it is proved in \28f that weak solutions are unique whenever H 
satisfies the structure conditions 


Hp{t,x,p) -p > rH{t,x,p) - 7 
\Hp{t,x,p)\ < /3(l-k|p|’'"^) 
H{t,x,p) > a\pY - 7 


for some r £ ( 1 , 2 ] and some positive constant a, fi, 7 . 

Therefore, in this case the convergence established in the above theorem holds 
for the whole sequence, and not only for a subsequence. 


4. A priori estimates. 

4.1. Norms and semi-norms. It is useful to define the following norms and 
semi-norms: 


^il carattere rispetto a x dove viene realmente usato ? Idem per la g 
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for a grid function v = we define 


lkllL'=(T2) 


lklllVC«(T2) 


^ ij 

+ {D+v,,,) 


2 \ 2 


Il*(T 2) + plvKi'*(T2)j ) 


where v and v are defined in (12.11) . We shall also write |r'|//i(T 2 ) 
ll^llffi(T2) = lkllvFi. 2 (T 2 ), and define the discrete scalar product: 


'i'J 


(4.1) 


(4.2) 

(4.3) 

kllVl>2(T2), 


We recall the discrete Sobolev inequality: for any s < oo, there exists a constant 
C such that for any grid function v, 

||?^||lo(t2) < C (^||r'||L2(T2) + |r'|_H-i(T2)) ■ 

For s > 1, we define the dual norm ||w|| w-i,s'c]r 2 )) 7 + f ^ 


,, ,, (^',^^)L2(T2) 

lklllV-i,*'(T2) — sup 

u)#0 ll^i’llM/i.s(T2) 

Define Qh,/\t = At{0,..., Nt — 1} x T^. For a function w defined on Qh,At: ^ = 
0 < n < Nt, we define for s £ [1, + 00 ), 


Nt 




(4.4) 


n—0 


4.2. First estimates. Hereafter, the constants appearing in the a priori esti¬ 
mates, for example c, C, are independent of h and At. In this paragraph, we state 
the first a priori estimates stemming from the structure of the system. Although we 
have already given the set of running assumptions, we think that it may be useful 
to specify which assumptions are really required by each particular result. 

Lemma 4.1. Under Assumptions (gij and (gs), if F is bounded from below by 
a constant F, uq is continuous on T^, then for all i,j,n, 


> u — T [ F — max H(x,0) 

' ~ xGT2 


(4.5) 


where u = min 2 ,gT 2 uo(x). 

Lemma 4.2. Under Assumptions (gi), (gs) and (gs), if F is bounded from 
below by F, uq is continuous on and mx is bounded from above by friT, then 
there exists a constant C such that 


Nt — 1 2 

9 q(.XiJ,[yhU^+^]^^.) <c, 

k—0 i,j 

N'T — 1 

h^Nt ^ ^5(a:.,„[V^u'=+i]^p<C, 

k—Q ij 

Nt — 1 

h^At 

fc=0 i,j 
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(4.6) 

(4.7) 

(4.8) 





Proof. Consider = nAtF{mT) for all i,j,n. We get immediately 


= F{mT) 


(4.9) 


Subtract (gH) from (1^ and multiply the resulting equation by ni^j — rriT- Sim¬ 
ilarly, multiply (1^ by ~ Adding the two resulting identities and 

summing with respect to n, one gets: 

Nx — 1 

{dqixi.j, .p • - g{xi,j, [V,iu'=+i] 

k—0 i,j 
Nt — 1 

+ h^At [V/.m''+^] . p 

k—0 i,j 
Nt — 1 

-f h^At ~ '^T){F{mlj) - F{fhT)) 

k—0 i,j 

=(m^^ - fhT,u^'^ - T F(TOT))L2(Tr2) - (to° - mr, u°)i2(Tr2). 

(4.10) 

1. Since — ffiT is nonpositive with a bounded mass, and since u" is 
bounded from below, see (|4.5I) . the term (m^'^ — rriT^u^'^ —T F{mT)) 

in the right hand side of (14.101) is bounded from above by a constant inde¬ 
pendent of h and At. 

2. It is straightforward to see that (to° — fhT,vP)L'^(j‘j^) < (1 + wt)||uo||oo- 

3. Since F is continuous, there exists a constant c such that F{f) < 
c, Vt > 0. Hence, 


N'T — 1 


N'T — 1 


k—0 i,j k—0 i,j 

4. Finally, ~ FiT)F{fhT) = T{1 - mT)F{rriT) 

From these observations, (14.61) . (14.71) and (14.81) follow from (14.101) and (13.IF □ 

4.3. A priori estimates from the discrete Fokker-Planck equation. 

The following estimates for the Fokker-Planck equation may have their own interest: 

Lemma 4.3. Assume (gi) and (gs). Let ip be a non decreasing and concave 
function defined on K+. For any grid function v = (vij), any positive grid function 
m = (mij) and any positive number g, 


^7)j (n,m)^/>(mij) < | ^[V?,?/'(w)]m ' hm]i,j 




1 f ^ '\ ^ 






(4.11) 







In particular, if m does not vanish, then for ipiz) = ln(z), 


'^■‘3 ^ 12) 

^,3 ' id 


Proof. By the definition of T, we can split the sum S = j %,j{v, m)'ip{mij) 
as follows: 

S = - '^m^jVqgixij, = Si + S2 + S'4, 




where 


1 ^9 


^i3 


1 ^9 

id ^ 

1 ^9 

S3 = --^'^mij-^{xij,[Vhv]ij){'tp{m^j+i) - ipiniij)), 

id ^ 

1 ^9 

It is enough to focus on S'! since the same arguments can be used for the other 
sums. Since g is nonincreasing w.r.t. < 71 , 

1 ^9 

Si < 'fl/hv\ij){i){mi+i^j) - V’(TOij))+. 

i.i ^ 

Since fi is nondecreasing, if mt+ij > uiij, the factor {'ip{mi+ij) — 'if>{mij))+ can 
be rewritten 


^z+l,j ^^id 

Since ip is nondecreasing and concave, mj+ij- > mij implies that 


0 < 




mi+ij) - ip{mij)){m^+ij - rrnj 

rriij impliei 


Hence, if mi+ij > mij, then 

{'lp{m^+lJ) - 'ip{mij))+ < ('ip'{mij){ip{mi+ij) - 'if>{mij)){mi+ij - rrnj, j , 

which implies that 


*^1 S Ip 








Lemma 4.4. Assume (gi) and (gs). Let tp be a non decreasing and concave 
function defined on 1R+. For any positive grid functions fc = 0,..., Nt, 

N'T N'T 

fe=n+l id k—n+1 id 

(4.13) 
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If G K-h for all k G {0,..., JVt} and does not vanish, then 
Nt — 1 

(4.14) 


k—n 




^,3 


Proof. Since ip is non decreasing, 

Nt Nt 


k—n-\-l i,j 


fe=n+l i,j 
^k-1 


From the concavity of if < m^j , then j) —— 
m'f~^). Then (14.131) follows from the last two points. 

Let us turn to (14.141) : for any e > 0, 

N'X' — 1 

X X Hmlj + e) =h? ^ m^J + e) “ X + e) 

k—n i,j 


^,3 




Nt 


- X X (Hmlj + e) - ln(mf_j- ^ + e)), 

k—n-\-l i,j 


and (14.131) with ip{z) = ln(z + e) yields 


Nt 

X 

k—n-\-l i,j 


Nt 


-h^ X + + X Xri^K7'-0+ 


, mj" ■ + 

k—n-\-l i,j ^i3 
Nt 

- X X("^h^ “ "^t)+ ^ 1’ 

k—n-\-l i,j 

where the last estimate comes from the fact that the grid functions all belong 
to ICh- Hence, 

Nt — 1 

X X ~ ^h) Hmlj+e) < X +e)-^^ X Hmlj+e)+l. 

k—n i,j 


^,3 


^,3 


and (|4.14|) is obtained by letting e tend to 0. □ 

Lemma 4.5. If G JCt and fgij (gs) hold, then there exists a constant C 
such that, for any number rj, 0 < rj < i/, a solution (to"^) of i2.9\) satisfies 


Nt — 1 




< C + h2^m^-|ln(mf7)| 


HHTl) 


maxh^ m", I ln(m’),-)| + (t^ — r])At |v7^ 

k=0 

U2 A -i ~^ 


k—0 ij 

For all a G [1,2)^ there exists a constant c such that 
ll’™llL“(Qh.At) 


(4.15) 


Nt — 1 

<c ( l + /i2X™f7’|ln(mf/)| + h2At X™bi 9qix^,J, 

k—0 i,j 


(4.16) 


Proof. 


10 

















Step 1. Take e > 0 and consider m". = to". + e. Note that to”. > 0 for all 
i,j,n. Multiply the second equation of (12.111) by ln(TO"j) and sum for all i,j and 
k = n,, Nt — 1: 


N'T’ — 1 


lyh'^At 


0 E E K. 1 ' - - 

k—n i,j 

N'j' — 1 

+ h'^At ^ y^ Ti,j{u''~'~^,m'") ln(TOf ,). 

k—n i,j 

From (I4.12|) and (I4.14L we deduce that 


k—n i,j 


U2\f Nt-1 

fc=n i,j 

^l + ^'E'^fo ln(^fo ) + ^ E 

ij k—n i,j 


and since m’l, ln(mj^ .) > —e 


A^t —1 


2,j A;=0 i,j 

U2A.J. ^T-l 2 

- )i + E E^bi| 5 ' 9 (^*o’■ 

2,ji ' i,j 


Consider now the quantity Wrh'^\'jji^j 2 y i.e. 


(4.17) 


\^\ hhti ) = E - v^) +E(\/*t+i-v/^) • 

Since > 0, we can write ^y^rhEiy “ ■\/^^E) ^ 

where (Il(''TO'')ij is defined in (12.IF Since the inequality (i/l + z— 1)^ < zln(l + z) 
holds for any number z > —1, we infer that 

- \f^) < h{D+ m%j ln(l + h ) 

'^•<3 

= h{D+7h^)ij (ln(TOf+i_^) - ln(TO,^^j)) 

= h^(i:)+TO'')ij(i:>)'' in(TO'=))ij. 

Since the same kind of estimate holds for ~ obtain that 

|\/^|^i(t2) < E ’ (4.18) 

id 

and the fact that to satisfies (14.151) follows from (14.171) and (I4.18F 
Let us now prove (I4.16F consider a G [1,2): there exists a unique number p 
1 such that ^ f we have the interpolation inequality ||TO^||ia(x 2 ) 

l|w'=||2i(T2jl|w''llip(T2)- But ||TO'=||ii(T 2 ) = 1 + e and = ||\/rhfc||i 2 p(T 2 ). 

From the discrete Sobolev inequalities, we deduce that 

||m'=||2c(T2) < (l + e)^ll^/j^ll22i>(T2) < C'll^^llffi(T2) < C'(l + l|\/^llffi(T=j) > 
which yields that to satisfies ()4.16p by summing for all k and using (j4.15F 
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lA IV 











































Step 2. We obtain that m satisfies (14.151) and (I4.16P by letting e tend to 0. □ 
Corollary 4.6. With the same assumptions as in Lemma for any a G 
[1,4/3), there exists a constant c such that 




m ^ — m 




At 


w-ro,(ji)) 


Nt — 1 


k—0 i,j 


^,3 


(4.19) 


Proof. Take a G [1,4/3). We start by observing that 

Nt — 1 

\\Dum\\2^^Q^^^^<Ch^At Y, Y\^tm% + \D+m%. 

k—0 i,j 

Let us estimate Si j by some quantity depending on 

Si ^ Holder inequality yields that 


N'T — 1 

EEi^ 

k—0 i,j 


-\-^k I cx 


'Nt-1 

^1 EEi^ 


k—0 i,j 


' N'T — 1 


+ . 2 (Df 


' Nt — 1 

EE 

ij 

' Nt — 1 


(L»+m'=^- ■ ^ 


)'t'^3 


= E Y^iD+m^)^J{D+h^{rn’^))^,J\ E E 


k—0 ij 

Standard calculus yields that 
therefore 


k—0 i,j 


{D{ 

{D+ ln(m'=)),j 


1 -# 


(J+ln 7 rl^)L- - 


Nt — 1 


E E 


+ ^fe la 


m 


k—0 i,j 
^ Nt — 1 


' Nt — 1 


<h^At\ Y E(^^"^^)bjPi'ln(TO'=))ij Y E(”^t 


A;=0 i,j 


k—0 i,j 


( Nt-1 

k=0 i,j 

Note that 1 < 2 ^ < 2. The estimate on \\Dh'>n\\‘fcc(^Q^ follows by using the 

same argument for bounding J2kN^ Si j then (14.171) and (14.161) . 

From (1^ . we deduce that for all grid functions w, 


h"^ Y 


- K,3 

At 




=vh?YNhm'^)i^j ■ {Vhw)ij +h^Y'N:,i9q{xi,jA^hu'^^\j) ■ \^hw\ij. 
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Hence 


^ra+l _ 

At 


,■ 0 ; ) 


/L2(t2) 


< J2|m"|vyi,<»(T2)|w|T^i,„'(T2) 




and 


Nt — ^ h-L^ h 

m ^ — m 


At J2 


k=0 


At 


W-^.-{Tl) 

+ j h^At ^ 9 q (xi,j, [V/,m'=+i] . J 




A^t —1 


(4.20) 


k =0 i,j 

Note that < 2. From (I4.16L 


Imp c 

i=-“(Qh,At) 




Nt — 1 


<c ( l + P^m,^^|ln(m^/’)|+PAt ^ X] 9 q(.x^J,[y 

k —0 i,j 


^,3 


The desired estimate on At J2k=o ^ 


m —m 


At 


vy-i>“(T2) 


(4.21) 

follows from (I4.20I) - (I4.21I) 


and from the estimate on \\Dh'm\\1a(^Q^ D 

Collecting the above results together with Lemma [4.21 we obtain the following 
conclusion: 

Theorem 4.7. If F is continuous and bounded from below by a constant^, if 
(giA (gsA (SsJ hold, if uo is continuous, then there exists a constant C such that 
a solution (u, m) of i2.12\} satisfies and for all a G (1,4/3); 


/ Nt-1 


fe=0 


j.^fe+1 _ jjjk 


At 


< C. 




(4.22) 


5. L^-compactness results. In this section we prove the L^-compactness of 
DhU whenever the discrete heat equation has bounded data. More precisely, we 
assume that u = {u'^)n=o,...,NT satisfies 

= /A (5-1) 

for all 0 < i,j < and all n, 0 < n < Nt, where the data / = (/";) ^.nd the 
initial conditions = (m°;) are supposed to satisfy 

‘S'IP°IIli(T 2) + ll/llz,i(Qh,At) — ^ (^-2) 

for some c independent of h and At. In what follows, we reconstruct functions 
on Q from the grid functions u, and we prove the convergence of these functions 
as h and At tend to 0, at least for subsequences. Lemma [Q] below is concerned 
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with piecewise constant functions built using u. It is similar to results that can be 
found in Gallouet et al, see e.g. [unnii] in the context of finite volume methods. 
Lemma [Q deals with approximations of the gradient with respect to x. It seems 
new to the best of our knowledge and may have an independent interest. 

Lemma 5.1. Let Uh^At be the piecewise constant function which takes the value 
u^'j^ in (tn,tn+i) X (ih — h/2,ih + h/2) x {jh — h/2,jh + h/2). There exists a 
subsequence of h and At (not relabeled) and a function u such that Uh,At u in in 
L^{Q) for all (3 e [1,2). Moreover, u e L“(0, T; IL1’“(T2)) for any a S [1, |), and 
there exist a bounded Radon measure ft in Q and a bounded Radon measure fo in 
such that u is the unique solution of 

jdtU-vAu = f inQ, 

\uiO,-) = f° inT\ 


Proof. Using the bounds on the data, we may show with the same argument 
as in S OI that 


lL'5(Qh,At) 


e Nt 

At ^ \Dhi 


n—1 




< C 


for any (3 G [1,2) and a € [1,|). From this estimate and (15.11) . we deduce that 




At 


- ||w-i.c(t 2) is uniformly bounded. 


Recall that u^ At is the piecewise constant function which takes the value u'fj in 
(t„, t„+i) X {ih — h/2, ih + h/2) x {jh — h/2,jh + h/2). We can apply the discrete 
Aubin-Simon lemma in |14j (Theorem 3.1): up to the extraction of a subsequence, 
Uh.At converges to a function u in L^{Q), and in fact in Ll^{Q) for all (3 G [1,2). 
Moreover, it G L“(0, T; 1U^’“(T^)) for any a G [1, |). 

Let fh,At be the piecewise constant function on which takes the values in 
{tn, tn+i) X {ih — h/2, ih + h/2) x {jh — h/2,jh + h/2). Up to the extraction of a 
subsequence, fuAt converges in the weak-* topology to some bounded Radon mea¬ 
sure f on Q. Call u/ the piecewise constant function on which takes the values 
u/ j in {ih — h/2,ih + h/2) x {jh —h/2, jh + h/2). We may assume that u/ converges 
to a bounded measure f^ on T^. In particular, testing (15.111 with smooth functions 
and passing to the limit, this implies that it satisfies 

/ / u{—ipt — vAip)dxdt= / / ipdf+ ip{0)df^, 

Jo jt 2 Jo jt 2 

for every tp G C'^{Q) such that ip{T) = 0. Notice that u is the unique solution of 
the above weak formulation. □ 

We now define an approximation of Dit from the grid function u. For a real 
number z, let floor(z) be the largest integer that does not exceed z, ceil( 0 ) be 
the smallest integer that is not less than z. Let Duh.At be the piecewise constant 
function from Q to which takes the value 




floor(-|),ceil(^) 




ceil(^),floor(^) 


value 
'{Dt 
{Di 


ii ^n+l) ^ 

(*!.(* +1)1) X 


1)|)' 

. More 

explicitly. 

Duh 

,At takes 

the 



in 


1 ^n+l) 

X 

(ih, ih + 

1) 

X 

(jh,jh + 

1) 



^ in 

{^n- 

t ^n+1) 

X 

(ih, ih + 

1) 

X 


■h) 



-i) in 

{tn-, 


X 

(ih — L, i 

ih) 

X 


■h) 



^ in 

{in: 

, tn+l) 

X 

(ih — I, i 

ih) 

X 

(jh,jh + 

!)• 
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Lemma 5.2. Up to the extraction of a subsequence, the functions Duh,At con¬ 
verge a.e. to Du in Q, and in L'^{Q) for any a € [1, |). 

Proof Since u is the unique weak solution of (15.3p , for every sequence of smooth 
functions fls and smooth initial data converging to fl and to ffi respectively, in the 
weak-* sense of measures, (which, for instance, can be constructed by convolution), 
the smooth solutions Us satisfying 

j dtUs - vMJs = fis in Q, 

\us{0,-)=il°s 

will converge to u, e.g. in L“(0,T; fL^’“(T^)) for any a G [1, |), see e.g. [7]. 

We now consider the finite difference approximation 

rrn-l-l _ jjn 

(5.4) 

with {U^)i,j = M°(a:ij) and = fi5{tn,Xij). 

Let (T be a positive real number: let be the piecewise linear function defined on 
M by Ta{z) = max(—cr,min(z,( t)). Let the grid function e on Qh,At be given by 
= uf j - Ug i j. Define also, for 0 < i,j < Nh, 

I” = {* s.t. maxde'VijI, KJ] < cr} ; = {j s.t. max[|e"j-+i|, |e"j-|] < a] . 

For any s : 0 < s < 1 we have 


Nt —1 Nt —1 

n—0 i,j n—0 i,j 


The first term can be estimated as follows 


Nt —1 Nt —1 

h^At ^ ^|D+e))+Y<^^At E E E iDtel-^^DtT,,{e-+\,r. 

n=0 i,j n=0 j igx"+'- 


Nt — 1 


+ l|s 


E E E \Dt<, 

"=0 J 


Nt — 1 

<Ti-f I h^At J2Dte:+^DtT,,{e-+\j 

n=0 i,j 


A^t —1 

+1 ^'At ^ Ei^^^n 

n=0 i,j 


Nt — 1 


^ Nt-1 

E E E 1 


1-s 


<T^-i I h^At Y, 

n=0 i,j 


Nt — 1 


+ u^At ^ Y\D 


-I- ra+li 
1 ^z,j \ 


'e7^lj\ + KV\ 


1-S 


E E 

n=0 i,j 


n=0 i,j 

Similarly we estimate the term with Di^ using the set J7)", and overall we deduce 
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that 


N'j' — 1 / Nt — 1 

h^At E \Dhel+^\^ < cTi-f h^At ^ ' DhTAe^+\, 

n=0 i,j \ n—0 i,j 

(5,5) 

for some constant c only depending on s. We estimate the first term from the 
discrete equation 



Nt — 1 iV'j’ — 1 


n—0 i,j 

Nt — 1 

-h^At ^ 

n—0 i^j 


n—0 i,j 


e"+l _ . 

’ 1 _ ALt 

At ^ 


which implies, using that {x — y)T„(x) > Q^ix) — Qa-iy) for the nonnegative and 
convex function 0(s) = T^{r)dr, 


N'T — 1 

n—0 i,j 

Therefore, we deduce from (15.5|) 


A^t —1 s 

n=0 i,j 

Taking the minimum of the right hand side w.r.t. ct, and using the bounds for 
fjLS, jjLos and the data in (15.2p . we see that 

for some c and 9 depending on s but not on h or 5. Recalling the definition of e, 
and the estimate on the discrete gradient, we have proved that 

II \Duu - < c||u - ■ 


Hence, 


II \DnU - < c||u,,At - Us,h.AtfL{\~Qt^,y 

where Uh.At has been defined in Lemma o and Us,h.At is the piecewise constant 
function that takes the value in (t„, tn+i) x (ih—h/2, ih+h/2) x {jh—h/2,jh+ 

V2)- _ ” _ 

Let us also define DUs.h,At from the grid function in a similar way as Duh^At- 
it takes the values 


{{DtU, 


+i'\ 

<5 ffloor(7 


),ceil(i) 


{D 


+ yn+l) 


ceil(^),floor(^) j 
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in X (i|, {i + 1)|) X (j + l)f)- Therefore, we see that 

II \DUh,At - £>m|'’||l1(Q) 

<11 \Duh,At - DUs^hAt\’'\\L'^{Q) + II \DUs,h,At - DUs\‘'\\l'^{q) + II \DUs - T)un|ii(Q) 
<c\\uh,At - ^S,h,At\fLi{Ql^^) + II \DUs,h,At - DUs\''\\l^Q) + II \DUs - Du\'‘\\l^q-^ 
<c (\\uh,At - '“llLi(Q;.^_it) + 11^ “ + ll^'5 - ^<5.?>,Ai|lLi(Qfe_^j)) 

+ II \ms,hAt - dusI^Wlhq-j + II \DUs - Dui^WmQ) 

where we have used (lOI) to obtain the third line. At fixed S, since Us is a smooth 
solution of the heat equation, the discrete approximation Us,h,At converges to Us 
in L‘^{Q) and DUs,h,At converges to DUs in L^(Q;R^). Using also Lemma lOl we 
get that 

hmsup|||Du/,,At - i:iu|'*||Li(Q) < c\\u-Us\f^lrQ) + \\\DUs - L>u|'*||ii(Q) . 

h,At^0 

We conclude using the strong convergence of IJs to u in L“(0, T; tU^’“(T^)) for any 
a G [1, |) (see e.g. [7], as <5 —>■ 0. So 

II \Duh,At-Du\^\\mQ)'^''\~"°0 

which in particular implies that Duh,At converges to Du a.e. in Q and then, by 
Vitali’s theorem, in L°‘{Q) for any a G [1, |). □ 

Remark 4. As a consequence of Lemma 1 5. ill for any ^ G K^f/iere exists a 
subsequence ofh and At (not relabeled) such that the maps (f, x) i—?> Duh,Atit, x+hf) 
also converge to Du a.e. and in L°‘{Q) for any a G [1, |). 

Remark 5. Alternative strategies can be used to construct a function defined 
on Q from the grid function u. For example, we can define Wh,At as the con¬ 
tinuous and piecewise trilinear function on Q which takes the values 
itn,Xij) and which is trilinear in the rectangles of the time-space grid Qh.At- The 
advantage of taking Wh,At instead of Uh,At is that the former has weakly integrable 
partial derivatives with respect to the spatial variable. Therefore, we can use directly 
Dwh,At instead of having to define an independent approximation of Du such as 
Duh,At- It is then possible to prove the following lemma, which may replace both 
Lemmas \ 5. 1\ a,nd 1 5. SI - 

Lemma 5.3. There exists a subsequence of h and At (not relabeled) and a 
function u such that Wh,At u in L“(0, T; 1U^’“(T^)) for any a G [1, |) and in 
T^{Q) for all jd G [1,2). In particular, Wh,At u and Dwh,At —>■ Du in L^{Q) and 
almost everywhere in Q. 

Proof. The strategy of proof is similar except that we may directly use the 
continuous version of the compactness lemma of Aubin-Simon, see j29j . for the 
function Wh,At- D 

6. Prom the discrete to the continuous system. 

6.1. A priori estimates and compactness. Let Uh,At and Duh,At be the 
piecewise constant functions defined in Lemmas 15.11 and 15.21 respectively: up to 
the extraction of a subsequence, we can assume that Uh,At u in L^{Q) for all 
(3 G [1, 2) and that Duh^At a.e. to Du in Q and in L°‘{Q) for any a G [1, |). 

Let mu.At be the piecewise constant function which takes the value mf ^ in (t„, t„+i) x 

{ih — h/2, ih + h/2) x (jh — h/2,jh + h/2), and Dmu^At be the piecewise constant 
function from Q to which takes the value 

{{D))m )flooqi)_(.eil(i) > )ceil(l),floor(i)) 
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in X + X + From Theorem 14.71 we may also 

assume that mh^At —>■ m in L^{Q) and almost everywhere in Q, and that Dmh,At 
Dm weakly in L°‘{Q) for any a G [1, |). Moreover, for all rj > 0, there exists a 
constant such that for all 2 > 0, + c^. This fact and estimate (14.81) 

yield the equi-integrability of F{mh At)- By Vitali’s theorem, F{mh At) —t F{rh) in 
L\Q). 

From the observations above, the piecewise constant function which takes the value 

_ yn 

- F{ml^) 

in [tn, tn+i) X {ih—hl2, ih+h/2) x {jh—h/2, jh+h/2) converges to ^ — vAu—F{m) 
in the sense of distributions. 

6.2. Stability of the discrete Bellman equation. We now pass to the limit 
in the discrete Bellman equation. 

The main difficulty is to handle the nonlinear term g{xij, here 

we wish to use the a.e. convergence of the gradients obtained in § O We adapt 
the method used for continuous problems in [8]. Note that ^ is the value 

taken by the piecewise constant function with values in 

( el • Duh,At{- + |el), el • Duh,At{- - fel), el • Duh,At{- + |el), el • Duh,At{- - |el) ) 

at {t,x) such that \xi — ih\ < h/2, \x 2 — jh\ < h/2, tn < t < t„+i. From the 
continuity of g, the consistency assumption and Remark HI 

9h,At -t g[x, Diu, Diu, D 2 U, D 2 U) = F[{x, Du) a.e. in Q, 

where gh,At is the piecewise constant function which take the value g{xij, . ^) 

for {t,x) such that \xi — ih\ < h/2, \x 2 — jh\ < h/2, < t < tn+i- 

Let now <^5 be a smooth function on such that > 0, with ip{T) = 0. We 
multiply the discrete Bellman equation by (/?(<„+!,and sum for all i,j and 
n = 0,..., Nt — 1- Since, by convexity, 


g{x,q) > g{x,0) + gq{x,0) ■ q 

the regularity of g w.r.t. x and the L^-compactness of Duh,At allow us to apply 
Fatou’s lemma obtaining 

Nt — 1 „ 

liminf^ ''^g{xij,[VhU^^^]. .)ip(tn+i,Xij) > / FI(x, Du)(p dxdt. 

n=Q i,j ’ 


Passing to the limit in the other terms of the equation, we deduce that 

— / uifitdxdt+u / DuDipdxdt+ H{x,Du)(pdxdt < / uo(f{0)dx+ F{rh)ipdxdt. 
■IQ -IQ -IQ dT2 Jq 

We now wish to obtain the reverse inequality, which is the difficult part. We start 
by noticing that, since the monotonicity assumption implies 

g{x,qi,q2,q3,qi) < 9{x,-q//,q^,-q^,qt) 


from (EU and (USD and the fact that g{x, 0) is bounded, we know there exists A > 0 
such that 


g(x, 91 , 92 , 93 ,94) < [1 + ( 9 i )^ + (9^)^ + (93 + (9^)^] . 
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( 6 . 1 ) 




We multiply the discrete Bellman equation by e ip[tn+i,Xij) and sum for all 

i,j and n = 0,..., Nt — 1- We obtain 


hr /St 


N'j' — 1 


n—0 


'i'J 


n+l _^n 

At 




N'j' — 1 

+ vh?At ^ 

n—0 i,j 
Nt — 1 

+ h^At 

n=0 i,j 


( 6 . 2 ) 


N'T — 1 

= h^At ^ 

n—0 iyj 


Since u is uniformly bounded below, the last term converges by dominated conver¬ 
gence, so 


Nt — 1 


ly-T — 1. p 

lim y (f(tn+i,Xij) = / F{fh)e~^'^ (pdxdt. (6.3) 

^ - rt ^ r J Q 


n—0 i,j ^ ^ 

By convexity of s i—>■ and since (/>(T, •) = 0, we have 


ii”t‘ - ti". 


E E ai 

n=0 i,j 


1 ^ e-^<J - e-^<V 

<-h'^At y y- — -(/7(t„+i,x,j) 

n—0 i,j 

= \ h^At y ' y ^itn+UX.,,N + 1 ^2 ^ g-A<, 

n—0 i,j 


A ^ ^ At 

n=0 i,j 

and so, again by dominated convergence, 

Nt-I n 

limsup h"^ At EE^ 




h^O 


Nt-I 

^'’,7 ^^5? —Au /, \ 

e Xjj') 


n=0 i,j 


At 


-Aiy 


dxdt + \ [ e (^(0) dx . 

A Jt2 


(6.4) 


Let us deal now jointly with the second and third term in (16.21) . First we split the 
energy term according to the sign of {u)ij (and respectively); indeed, 

we can write 




g _ g 


— Xu^‘ 


") - iDt{u-+NN ( 


—/^—\ur 


+i,i — e 


— Au” 


zzyu"+i)..,i^+(e- 
= (Z^+(u"+i)..,)+ ( 

+ (z?y«”+^)..,)+(^'vy - - (z?y«”+')..,)-«+y - 

and the same for the term with D 2 ■ Reordering the indexes in the sum, this means 


■)< 


n-\-l 

3 
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that the Di part in the second order term can be read as 


A^t —1 




->'<T - \ 


n—0 i,j 
Nt — 1 




- e 


-A«"+l \ , n+1 




-lyh^At j2 ^ (e 

n=0 i,j 
Nt — 1 

+ ,yh^At J2 J2^D+{u'^+\,)+D+iip^+\,e-^<-" 

n—0 ifj 
N'j' — 1 

n—0 ij 


which is equal to 


N'T — 1 


i^h^At J2 

n—0 i,j 


_A«"+i -A«"+7 

t,3 — g t-1,. 


Nt — 1 


lyh'^At 


n+l _ n+1 

^i,j ^i-l,j 

— Xu^!^} . —Xu^ 

e •+n3 — e •■■ 


■(f. 


n+l 

'■ij 


,,"+1 _ „"+i 

“i+i.i ^i,j 


n+l 


n=0 i,j 

Njp — 1 


+ ^ h-^At X X ' - ++(r."+i).+-e-^<+c. } i?+(++i) 

n=0 i,j 

We proceed similarly for the part with D 2 . Therefore, 

N'J' — 1 

vh^At XX • Dh{e-^<T j{tn+ux,,jN, 

n=0 i,j 

N'J — 1 

+ h'^At Y, X:s(x.„[v.„"«].ye--SV(i„+„.,j) 

n=0 i,j 

Nj — 1 


= ixh^At X X - {D+iu^+N,)-e-^<+^"-}D+{j^+Nj 

n—0 i,j 
Nj — 1 

+ ixh^At X 


n=0 i,j 

Nt-1 e->'KV - ^->'K-ii 

' h'^At y : y ; i ( D +(„”+')._ i ,,)+'“^^^^—^^,.-+1 


n=0 i,j 
Nj — 1 


,," + 1 _ „" + l 
-A«"+b. „-A«"ti 


'h^At X „rr'_Er^ 

'^2+1,j 

_An^+i -An^t^ 


n=0 
iVx- — 1 


'hAt X ' h.i-ll ' I n+l „n+l 


,,n+l\ . 'i+|2' 


n+l 


n—0 i,j 
Nj — 1 


->'<V 

^.J + l — g 


• X J2\(^tiu"^Njyr—j+i ——+, 


n+l 


n=0 i,j 
Nj — 1 

+ h‘^At X Y 9(^+3 

n—0 i,j 
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The first two terms in the right-hand side converge to v fj' e ■ d(j) by 

Lebesgue theorem, since Duh,At converges strongly in L^, (p is smooth and 
is uniformly bounded and a.e. convergent. As far as the remaining terms are 
concerned, using that 

, ,, 

_ _ _ _ ^ —Am.ax(s,s ) 

s — s' ~ ’ 

and due to (I6.1I1 . we observe that the last five terms under summation are bounded 
above, so that we can again apply Fatou’s lemma, on account of the a.e. convergence 
of Uh,At and Duh,At- Therefore, we conclude that 


Nt — 1 

limsup vh'^At ^ ^ (p{tn+i,Xij))tj 

n=0 i,j 

N'j' — 1 

+ h^At 

n—0 i,j 

< i’ j DuDp dxdt — V \ / \Du\^ p dxdt + H{x,Du)e~^''p dxdt. 

JQ Jq 

(6.5) 

Putting together (I6.3I) - (I6.4I) - (I6.5I) . we deduce from (16.21) that u satisfies 

-- f e~pt dxdt + — ( p{{)) dxV f DuDp dxdt 

^ Jq ^ Jt^ Jq 

— vxf p dxdt + ( H{x, Du)e~^''p dxdt > f F{m)e~^''p dxdt 

Jq Jq Jq 

( 6 . 6 ) 


for every smooth p > 0. In order to conclude, we need now to get rid of the 
exponential in the above inequality (16.61) . To this purpose, we first observe that 

G L'^{0,T;H^{T^))nL°°{Q). (6.7) 

This can be easily proved obtaining an a priori estimate on Indeed, 

whenever u is a grid function which solves for some data satisfying (lOD , we 
have 


h^At 


A^t —1 

EE 

n—O i,j 


At 


Nt — ^ 


^«+i) +,.h^At E 


n—O i,j 


<c\m 


for any bounded real function ^p(r). In particular, if ^l} is nondecreasing, this implies 

N'J' — 1 N'J' — 1 

n—O i,j n—O i,j 

where 'P(s) = tp{r)dr. Thus, since |'If(s)| < ||'!/'||oo|s|, one gets 

N'J' — 1 

i^h^At ^ ^i4,<+i.D,^«+i)<CMoo 

n—O i,j 

where C only depends on the L^-norm of the data. This is the desired a priori 
estimate; from which, using Fatou’s lemma, we deduce 


IQ 


\Du\'^ ip'{u) dxdt < C||^/>||c 
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On account of the fact that u is bounded below, we can take for example ^/’(r) = 
1 — to deduce that G L^(0, T; iJ^(T^)) for any /r > 0. 

Thanks to (I6.7|) , inequality (16.61) holds true not only for smooth functions ip but 
also for ip G H^{Q) fl L°°, through a standard density argument. Moreover, there 
is no loss of generality in assuming that uq G i7^(T^), so we extend u for negative 
t as identically equal to uq. Then, we choose 

ip{x,t)=^{t)y [ 6^^'= ds 

h Jt-h 


where Tfe(r) = min(r, fc) and ^ G Cl[0,T). Using the monotone character of s i—>■ 
we have (see Lemma 2.3 in i) 


lim sup 
/t-s-O 

< - 


e dxdt + 


A 


it dr 

IQ Jo 



V3(0) dx 





Moreover, ^\Tk(u)(x,s) converges to in L^(0, T; iJ^(T^)) and weak—* 

in L°°{Q), so we can pass to the limit as h —>■ 0 in the remaining terms of (16.61) . 
Finally, we obtain 


p pU p pUQ 

- it / ^(0) / dr 

Jq Jo Jt^ Jo 


+ u dxdt-uX 

Jq 


{u>k} 


+ [ idxdt> [ f^xdt. 

Jq Jq 


We conclude by letting k —>■ oo, thanks to the dominated convergence theorem: 


— / u£fdxdt- 








JQ Jq Jq Jq 

for every ^ > 0. Since the reverse inequality was already obtained previously, in the 
end, we proved that u solves 



uoi{0)dx+v / DuD^dxdt+ H{x,Du)^dxdt 

Jq Jq 


F{rh) i dxdt 


for every ^ G C',(([0,r)),^ > 0, and therefore for every This concludes the proof 
that M is a weak solution to the limit equation. 


6.3. Stability of the discrete Fokker-Planck equation. We now pass to 

the limit in the discrete Fokker-Planck equation. _ 

By (14.6L the L^-compactness of mh,At and of Duh,At, we deduce the strong 
convergence in L^{Q) for the piecewise constant function which takes the value 
m2jVqg{xij,[\7hu'^~^^]i,j) for {t,x) such that |a;i — ih\ < h/2, \x 2 — jh\ < h/2, 
tn < tn+i- Moreover, by the consistency assumption we have. 


Nt — 1 


n=0 ^^j 


/ Vqg{x, Diu, Diu, D 2 U, D 2 U) ■ {Diip, Diip, D2ip, D2ip)dxdt 

Jq 

f dH 

= / m-^{x, Du) ■ Dip dxdt 

Jq 9P 
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Therefore, we can pass to the limit in the weak formulation and deduce that m is 
a weak solution of the Fokker-Planck equation. 

We notice that the reg\i\a.rity Du)Du —H{-, Du)] G L^{Q) follows from 

inequality (I4.10I1 . by using Fatou’s lemma. Moreover, we also find that Du)]"^ G 

L^{Q). The regularity u,m G C°([0, T]; L^(T^)) follows from properties of weak so¬ 
lutions, see [28] . 

Finally, this concludes the proof of Theorem 13.II 
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